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p [ Resonant cyclotron scattering of the surface radiation in the magnetospheres 

of neutron stars may considerably modify the emergent spectra and impede ef- 
forts to constraint neutron star properties. Resonant cyclotron scattering by a 



ABSTRACT 



c§ ! non-relativistic warm plasma in an inhomogeneous magnetic field has a num- 

ber of unusual characteristics: (i) in the limit of high resonant optical depth, 
^ ! the cyclotron resonant layer is half opaque, in sharp contrast to the case of non- 

resonant scattering, (ii) The transmitted flux is on average Compton up-scattered 
by ~ 1 + 2[3th, where [3th is the typical thermal velocity in units of the velocity 
of light; the reflected flux has on average the initial frequency, (iii) For both 
the transmitted and reflected fluxes the dispersion of intensity decreases with in- 
creasing optical depth, (iv) The emergent spectrum is appreciably non-Plankian 
while narrow spectral features produced at the surface may be erased. 

We derive semi-analytically modification of the surface Plankian emission 
due to multiple scattering between the resonant layers and apply the model to 
anomalous X-ray pulsar IE 1048.1-5937. Our simple model fits just as well as 
the "canonical" magnetar spectra model of a blackbody plus power-law. 
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1. Introduction 

The magnetospheres of neutron stars are filled with plasma, with the minimum magne- 
tospheric plasma density being on the order of the Goldreich- Julian density ncj = B-f2/27rec 
(B is magnetic field , Q is angular velocity, e is the elementary charge and c is the velocity 
of light), Goldreich & Julian (1968). The magnetospheric plasma may efficiently scatter 
and absorb surface emission through resonant cyclotron scattering. If the plasma density in 
the magnetosphere equals the minimum Goldreich-Julian density, then resonant scattering 
is not important in the optical- X-ray range (Mikhailovskii et al. 1982, Eq. (7)). On the 
other hand, under certain circumstances the real plasma density in the magnetosphere may 
exceed this minimum value. In rotationally powered pulsars pair excess on closed field lines 
can be due to pair creation by high energy photons {e.g. Wang et al. 1998). In magne- 
tars {e.g. Woods & Thompson 2004), large scale currents flowing in the magnetosphere 
result in particle densities much larger than the Goldreich-Julian density (Thompson et al. 
2002). In accreting sources {e.g. neutron stars in X-ray binaries or isolated neutron stars 
accreting from the interstellar medium, ISM) the external plasma may penetrate inside the 
magnetosphere due to the development of instabilities in the boundary layer (Ghosh et al. 
1977; Arons & Lea 1980). Finally, the presence of a dense magnetospheric plasma is well 
established in the case of the binary pulsar PSR J0737-3039B, where the multiplicity factor 
(the ratio of plasma density to Goldreich-Julian density) was inferred to be ~ 10^ (Lyutikov 
& Thompson 2005). 

The presence of a relatively tenuous plasma may lead to large magnetospheric optical 
depths for scattering and absorption since in soft X-rays the resonant cyclotron cross-section 
is five to six orders of magnitude larger than the Thomson cross section, Eq. (6). Due to the 
large resonant cross-section, the amount of electron-ion plasma that must be suspended in 
the magnetosphere to produce an optical depth of the order of unity is tiny by astrophysical 
standards, of the order of 10^ grams total, Eq. (9). Both the emission patterns and spectral 
properties of the surface radiation are then modified in the magnetosphere, as photons are 
scattered and their energies are Doppler-shifted in each scatter. 

In this paper we describe salient properties of resonant scattering in neutron star mag- 
netospheres and investigate semi-analytically the influence of scattering on emergent spectra 
in the non-relativistic limit, when the plasma temperatures are smaller than mc? jhs (where 
m is the electron mass and A;^ is Boltzmann's constant). We employ methods of radiative 
transfer that were developed for line-driven winds (Chandrasekhar 1945; Sobolev 1947) and 
were previously applied to cyclotron scattering by Zhelezniakov (1996). We adopt a simple 
one-dimension model for photon propagation (called the Schwarzschild- Schuster method) 
which allows a clear description of the underlying physics. More involved Monte-Carlo sim- 
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ulations which take into account the full three-dimensional structure of the magnetosphere 
as well as relativistic effects are underway. 

The majority of prior work on resonant cyclotron scattering considered either the atmo- 
sphere or the accretion column of neutron stars (Wasserman et al. 1980; Nagel 1981; Lamb 
et al. 1990; Dermer & Sturner 1991), in which case the thickness of the plasma column H 
(centimeters or meters) was assumed to be much smaller than the thickness of the resonant 
layer (see Section 2) Ar ~ f3thr (r is the distance from the center of the neutron star and j3th 
is the thermal velocity in terms of the velocity of light). In this case the transfer of resonant 
radiation occurs in a quasi-homogeneous magnetic field. Higher up in the magnetosphere H 
becomes larger than Ar, so that the effects of the magnetic field's inhomogeneity become 
important. In this case properties of resonant scattering are different from the case of a 
homogeneous magnetic field and are somewhat counterintuitive. For example, Zhelezniakov 
(1996) has pointed out that in a strongly inhomogeneous magnetic field in the high optical 
depth limit resonant layers are semi-opaque. Following up on this work, we investigate spec- 
tral properties of the transmitted radiation. We show that if a considerably dense, warm 
plasma is present in the magnetosphere (so that both the resonant cyclotron optical depth 
and the frequency Doppler shifts during scattering are no n- negligible) , then (i) the trans- 
mitted spectra are modified, so that thermal surface emission has a resulting non-Plankian 
spectrum] this effect is most pronounced at optical depths ~ few, while for very large opti- 
cal depths the spectrum remains nearly Plankian; (ii) transmitted radiation is preferentially 
Compton up-scattered in frequency, so that the real surface temperature may be smaller than 
the one inferred from spectral fits; (iii) surface spectral features may be erased; (iv) scatter- 
ing in the optical and X-ray bands occurs in different regimes due to different polarization 
properties of normal modes (circularly polarized in optical and linearly polarized in X-rays) . 

We assume that plasma is tenuous, so that typical frequencies are much larger than 
the plasma frequency, and the typical wavelength is much larger than the skin depth. The 
refractive indexes are then close to unity and Debye screening is not important in calculating 
the scattering probabilities. Collisional processes may also be neglected for the densities and 
temperatures involved. We also assume that initially radiation is not polarized; this is a 
considerable simplification since neutron star surface emission is expected to be polarized, 
especially at magnetar field strengths (Lai & Ho 2003; Ozel 2003). 
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2. Resonant Scatter in the Magnetosphere 

2.1. Resonant Layers 

If the magnetosphere of a neutron star is populated by warm (non-relativistic) plasma, 
quasi-thermal radiation coming from the surface of the neutron star may experience strong 
scattering at the cyclotron resonance, located for a given photon energy at a radius 

^~^„J^)"^8.„.6V3(^)-"^ (1) 
\ emc ) VlkeV/ 

where r^s ~ 10^ cm is the neutron star radius, B^^s is the surface magnetic field, h = 
Bjsis/ Bcr-i Bcr = 4 X 10^^ G is the critical magnetic field, h is Plank constant divided by 27r, 
m is the electron mass, e is the energy of a photon. The importance of resonant scattering 
may be characterized by the resonant optical depth (Mikhailovskii et al. 1982; Zhelezniakov 
1996) 
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where a' is the angle between the incoming photon and the local magnetic field, ujb = eB/mc 
is cyclotron frequency, ares is a resonant cross-section 

_ ctt (1 + cos^a')cj2 

~ T(a;-c^s)2 + rV4' ^ ^ 

r = 4e^u;^/3meC'^ is the natural width of the first cyclotron line, n is the plasma density 
and we assumed dipolar fields with \dlnu!B/dhir\ = 3. The meaning of optical depth tq 
is that only e~^° fraction of photons coming from direction a' passes through the resonant 
layer without scattering. This is different from the relative number of photons that passes 
through the resonant layer after experiencing multiple scattering in the layer. 

In the spectral range of interest absorbed photons are almost immediately re-emitted 
since the synchrotron transition times ~ l/F are shorter than the dynamical time r/c for 

/4R2 4 6 \ 1/5 

-■^(^^ ~7.6xl0V/=cm (4) 
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Only microwaves and longer waves suffer true absorption; all the higher energy radiation is 
scattered. 

The resonant optical depth (2) greatly exceeds the non-resonant optical depth for Thom- 
son scattering tt ~ Rnax, 



1 keV\ 



10^ — 10^ (6) 



~ ~ lU I — I 

Tt 8re(jJB \ nujB J 

where = /mc^ is the classical radius of an electron. Estimating the total Thus, the 
presence of even a tenuous plasma in the magnetosphere may provide a substantial optical 
depth to resonant scattering. The corresponding optical depth is 

so that minimal magnetospheric density is not sufficient to produce considerable absorption, 
except near the light cylinder. In order to produce an optical depth of the order of unity the 
required density is 

?>B 3b mcmc^ r%q 
n ~ ~ — (8) 

Estimating the total volume of the scattering material as F ~ Anr^/S, the total number of 
scatterers required is only 

N ^!i5<iQ33^ (9) 

If the plasma is electron-ion, then this corresponds to a mass of only ~ 2 x 10^6 g. 



2.2. Resonant Scattering in an Inhomogeneous Magnetic Field 

magnetic field ^ the spatial resonance width is 

~ n 1 T = ^thLB ~ /9rr/3 (10) 

where Lb = —1/drlnujB and the last equality assumes a dipole magnetic field. 

Depending on whether the thickness of the scattering layer H is larger or smaller than the 
resonant width Ar = Pth^B, the radiation transfer occurs in two difi^erent regimes: a weakly 



^The magnetic field can vary in direction and magnitude. For a non-relativistic plasma the variation 
in magnitude is more important than the variation in direction, except in a narrow region near a' n/2, 
(Zhelezniakov 1996) 
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inhomogeneous field, PthLs ^ H (this is relevant to radiation transfer in photosphere and 
polar column), and a strongly inhomogcncous field PthLs H, relevant for radiation transfer 
in the magnetosphere (Zhelezniakov 1996). In a weakly inhomogeneous field a resonant 
photon that has been scattered once and whose frequency has been modified due to the 
Doppler effect, has a final frequency that typically lies outside of the resonant layer, so that 
the photon is likely to leave the system. For a continuum spectrum falling on the resonant 
layer this leads to formation of a narrow absorption line. In a strongly inhomogeneous field, 
Ar <^ H, a. scattered photon typically remains in resonance and can experience multiple 
scattering. If a warm plasma, Pth < 1, hUs a considerable fraction of the magnetosphere, 
the resonant width is smaller that the thickness of the plasma column, if ~ r, by a factor 
of Pth/ 3, so that resonant scattering occurs in the limit of a strongly inhomogeneous field. 
For relativistically hot plasma, f3tfi ~ 1, the width of the resonance becomes comparable 
to the inhomogeneity scale and the approach developed in this paper is not valid anymore. 
Resonant scattering in the relativistic limit has been studied in the context of early models 
of GRBs and pulsar magnetospheres {e.g. Lamb et al. 1990; Brainerd 1992; Sturner et al. 
1995; Lyubarskii & Petrova 2000; Gonthier et al. 2000) 



2.3. Polarization in resonance 

In neutron star magnetospheres both plasma and quantum effects influence the prop- 
agation and wave absorption through polarization of normal modes. Depending on the 
polarization of normal modes, there are two possible regimes of cyclotron radiation trans- 
fer. For weak magnetic field and/or low frequencies, plasma effects dominate over quantum 
effects. In this case waves propagating along the magnetic field are circularly polarized (if 
plasma is electron- ion) and in the case of warm plasma, (5th ^ 1; radiation transfer occurs on 
the extraordinary wave, while cross-sections for scattering of the ordinary wave is smaller by 

Alternatively, in a strong magnetic field and/or for high frequencies, wave polarization 
is determined by quantum vacuum effects. In this case waves are linearly polarized (for 
parallel propagation) and both modes are strongly scattered. 

In resonance, waves remain linearly polarized (dominated by quantum contributions) if 
the following condition is satisfied (Zhelezniakov 1996) 



Assuming that the optical depth is of the order of unity, uj^ ~ ujbc/t, dipole geometry and 
(3th = 0.3 we can find the photon energy above which polarization is determined by quantum 
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vacuum corrections 



Thus, starting from far UV, quantum vacuum corrections to wave polarization dominate over 
plasma contributions, so that normal modes are linearly polarized and have similar resonant 
cross-sections (Thompson et al. 2002). 



3. Radiation transfer 

3.1. General relations and Schwarzschild-Schuster method 

The photon transfer equation is best derived using covariant formulation (Lindquist 
1966). The invariant transfer equation for photon occupation numbers n is 

k^d.n = J Sff^ (13) 

where 5" represents a collision rate between photons with occupation numbers n and electrons 
with a distribution function /. Integration is over electron momenta d^p, 7 is the electron 
Lorcntz factor, A;^ is the photon four-momentum. Wc wish to write down the equation of 
radiative transfer for resonant scattering taking into account only linear terms in (3. This is 
best done in the frame where the electron is at rest and then transforming to the fluid frame 
(Hsieh & Speigel 1976; Blandford & Payne 1982). 

The invariant collision integral is 

.(..i;.P) ^ .0 / (( -^-K-'o^-^-'o) ) (^)= (^) II - 

\ In^ ) ^(^o,lo)(l + n(i/o,lo) 1 (14) 

where prime denotes a scattered photon, and subscript denotes quantities measured in the 
electron's rest frame. 

We are interested in terms linear in the thermal velocity of the particles, and thus we 
neglect recoil effects. In this case v'q — uq and the scattering cross-section is symmetric 
between the two states: 

diiQ J \ d\ Iq J mc \ \ ) J 

(15) 
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where //q is the cosine of the angle between the photon momentum and the direction of the 
magnetic field. Neglecting induced process, n <S 1, we find 



i^B I fe — / dQ'Q{l + i^l){l + i^'f){n{i^Q,l'o)6{iJo-iJB) -n{i^o,h)S{iyo-J^B)) (16) 
mc J J 

Transforming to the lab frame, vq — ^{1 — and using 

n(i/o, lo) = n{v\ V) ^ n{v, V) + ((1' -\) ■ (3) i^dM^, 1') (17) 

we arrive at 

k^d.n^UB^ [ /e^(l+/X^) / dn'{l+f^'^) X 

mc J ^ J 
((1 + ((!' vd,) n{v', ij!)b{v'^{\ - Hp!) - vb) - n{v, ix)b{v^{\ - - vb)^ 

where ji^ and ji'^ now should be expressed in terms of laboratory frame quantities, fiQ ~ 
(/X — /3)/(l — /3/i) and similarly for //[,. Assuming a non-relativistic one-dimensional electron 
distribution, /g = /(p||)(^(p±)/27rp_L, 7 ~ 1, and a stationary plasma (so that odd power of 
(5 average to zero), we find 

mc J J 
{n{u', i^')5{u'^il - f3i,') - vb) - n{v, ^Ji)5{v^{l - fJ/,) - ub)) (19) 

In order to investigate the salient features of resonant radiation transfer, wc make 
a simplifying onc-dimcnsional approximation for photon propagation (this is called the 
Schwarzschild-Schuster method in radiation transfer). More advanced analytical methods, 
e.g. , using Gauss formula for numerical quadratures (Chandrasekhar 1944) lead to compli- 
cated equations. [Note that Chandrasekhar (1945) used Gauss method for similar problem 
of line-driven wind, but neglected aberration effects and kinematic enhancement of cross- 
section, both linear terms in particle velocity In the Schwarzschild-Schuster method the 
integral in the source function is replaced by a sum over forward and backward propagating 
photons. Eq. (19) is then replaced by a system of two equations for / = ±1: 
22 /■ 

Id.ni = Y-^B / /e(/3)o?/3 (^('^^ ^')'5(^'7(1 - /3/') - vb) - n{v, V)b{v^{\ - pi) - ub)) 

(20) 



3.2. Radiative transfer in an inhomogeneous magnetic fields 



In an inhomogeneous magnetic field radiation transfer occurs both in real space as pho- 
tons diffuse through the resonant layer, and in frequency space as their frequencies are 
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Doppler shifted at each scattering. Neglecting the natural line width F, at each scatter there 
is a conserved quantity, the velocity of the particle 

^^o^^cy-c^ (21) 
uj cos a uj' cos a' 

Thus, resonant radiation transfer on particles with different velocities occur independently 
(Zhelezniakov 1996). This allows a simplification of the problem if we consider (3 instead of 
uo as an independent variable. Accordingly, we transform to variables x — /? by multiplying 
Eq. (20) by — /3/x) — ub) and integrating over v: 

(a^ + dr In VBdjj) n+ = '^^f'^^f^) - n+ + (3{n- + n+)) 

2 2 

(-9, + dr In VBdp) n_ = |^/e(/3) {n- - n+ - /3(n_ + n+)) (22) 

where n± = J di/n (i/) 5 (i/(l ^ — ub). 

In a weakly inhomogeneous magnetic field we can approximate cub — <^Bo ~ x/ Lb)- 
After dimensionalizing, x — > x/Lb, system (22) becomes 

x ■ 

{dx + dp) n^ = (n_ - n+ + /3(n_ + n+)) 

(-6>, + dp) n_ = ^ (n_ - n+ - /?(n_ + n+)) (23) 

where r^es = LBfe{P)'n'^e^/mc. 

System (23) differs from the similar equations of Chandrasekhar (1945) and Zhelezniakov 
(1996) due to the presence of velocity terms on the right-hand side. These terms come from 
the kinematic factor 1 ± /3 in the collision rate. Since we are interested in the non-relativistic 
limit /? < 1, in what follows we neglect terms with (3 on the right hand side. This allows 
us to obtain an analytical solution for the transfer function. In appendix A we verify the 
validity of this approximation using Monte-Carlo simulations. In the limit /3 — > system 
(23) becomes 

{dx + dp) n+ = ^ (n_ - n+) 

(_a, + a;3)n- = ^(n_-n+) (24) 

Resonant optical depth Tres is generally a function of both x {e.g. through the depen- 
dence of the density on x) and (3. For simplicity we assume that the density is constant, 
'Tres — 't{(3). Introducing new functions / and g {cf. Chandrasekhar 1945) 

n+ = e- / ^_ ^ e- / (25) 
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Eq. (24) becomes 

(d. + dfs)f^^-fg, (d^-dp)g^-^-ff (26) 
Which may be expressed as a second order partial differential equation for / only 

{dl -d})f + dp In Tres {d. + d0)f + ^f^O (27) 

This equation describes the transfer of resonant cyclotron radiation in a hot plasma in an 
inhomogeneous magnetic field. 

As a simple analytically tractable case we assume that the distribution of particles is a 
"water-bag" 

rW^lf- ^7J'*I*«' (28) 
0, otherwise. 



Eq. (27) then becomes ^ ^ 



Assuming that ujq is the frequency resonant at a; = 0, a;o = ^Bq, the resonance is given 
by a; = ^^^(l — ^) and ui = LJBo{l + r])- where rj = P — x, ^ = x + (3 are forward and backward 
propagating characteristics. Changing to variables ^,r), Eq. (29) becomes 



-^d^dj + T^/ = (30) 
The Green's function for this equation is 

GiC -Co,V- Vo) = /o (^-^^V{C-Co){v-m)^ (31) 

The initial condition for this system is that a monochromatic wave falls onto the inner 
edge of the resonant layer and that there is no radiation falling on the outer edge. In 
coordinates {/?, x} the inner boundary of the resonant layer is located at /9 = —Pth, x = —Pth, 
where r/ = 0, ^ = —2Pth (see Fig. 1). At this point we have n+ = S{u!—u!bo) — ^{^Bo{x—P)) — 
5{r])/ujBo- And ri- — Q oa g ca d^f = 0. Thus, the two boundary conditions are 

/(e = -2Pth, V) = S{rj) = -'^Pth, V) = (32) 



^This is a wave equation with massive scalar field, or equivalently, a telegraphist equation. 

^Zhelezniakov (1996) have neglected the term inside the resonance which is incorrect since the width 
of the resonance in dimensionless a; is 2f3th, the same as in /3 
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where the intensity has been normahzed to unity. 

The solution of Eq. (30) given the boundary conditions (32) can be found using standard 
methods of mathematical physics {e.g. Courant & Hilbert 1953). 



^ = + ^ \/ (jh Vi^-^o){v-Vo)) (33) 
8pth y rj-Tjo \4:(Jth J 



and 



^ g-(r?+e+2/3tt)ro/(8/3ih) 



n- = ^e-(^+«+^^-)-/(«^-)/o (^V(e-Co)(r/-r?o)) (34) 

In the case of a water-bag distribution 770 = 0, ^0 = —'^Pth, and we find the transmitted and 
reflected fluxes 



= it'-"'" (ii^V(2ft^-«)« + 2«) (35) 

where the transmitted flux should be evaluated at ^ = '^Pth ~ while the reflected flux 
n_ at 77 = 2Pth — ^- Returning to frequencies lu, 

/'/o {2^^ V^o(l + 2A^) -u;){u;- u;o{l - 2A/.))) (36) 



These relations are the sohitions for the transfer function, giving the fluxes of transmitted 
and reflected radiation for an initially monochromatic wave (see Fig. 2). 

The reflection probability, found by integrating (36) over frequencies, is 

1 - e-^ , , 

Prefl = ^ (37) 

In the large optical depth limit Prefi = 1/2. (Zhelezniakov 1996). 
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For small optical depths 



To 



To 



4A 



th 



(38) 



while in the limit of large optical depths 



^ 1/4 / 



l + A(3th 1 



-3/4 



(a;-a;o{l+2/3th))2ro/(16/3,\a;g) 

^ fl + - ^)(^ ^ (1 - 2/3,,)))) e(--o)^-o/(^«/^*>o) 

4v27ra;o V t^o t^o / 



(39) 



3.3. Qualitative properties of resonant transfer function 

The properties of the transmitted and reflected fluxes (36) are very different from Thom- 
son scattering and from resonant cyclotron scattering in a homogeneous magnetic field . 
(Fig. 2). In the small optical depth limit, tq <^ 1, the reflected wave has an approximately 
constant spectral flux, oc tq, centered at ujq with dispersion Acu/cu ~ 2Pth- The transmitted 
radiation, oc Tq , is up-scattered, spanning cuq < cu < {1 + APtfi)^o- 

In the large optical depth limit tq » 1 the spectral fluxes of the reflected and the 
transmitted radiation are equal, so that a resonant layer with high optical depth is half 
opaque, Eq. (37), see also Zhelezniakov (1996). This is in sharp contrast to the case of 
non-resonant scattering and to the case of resonant cyclotron scattering in a homogeneous 
magnetic field, when the transmitted wave is only oc e~'^ <^ 1. Half opaque resonant layers 
are also reminiscent of the case of line-driven winds (Sobolev 1947). The transmitted radiation 
is Compton up-scattered by a factor ~ 1 + 2(3th, while the reflected waves on average have 
UJ ~ u!o. Note that this up-scattering is done by a non-relativistic, static warm medium in 
the limit where recoil effects are neglected, (c/. Kompaneets 1956). Next, as optical depth 
increases, the transmitted and the reflected fluxes have increasingly narrower distributions. 
The typical dispersion around the central frequency is Au/ujq ~ 2-\/2/3t/i/-\/T for r^es ^ 1- 
This again can be contrasted with non-resonant Thomson scattering, in which case dispersion 
increases with increasing r^es- 

Consideration of photon residence time inside the resonance layer and number of scat- 
tering in the limit of high optical depth also bears a lot of surprises. Residence time inside 
the layer is independent of optical depth, while number of scattering scales as ~ r^es! To 
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understand this, consider a photon moving outward through a resonant layer. Suppose it is 

scattered backward at a point —Pt < x < [3t on an electron with resonant velocity (5 — x 
(in dimensionless units). It has already traveled a distance (St + x inside a resonance layer. 
After the scattering the photon energy is oc 1 + 2x and the new resonant surface is located 
at 2x, so that during scattering the photon is located a distance jSx — x from the edge of 
the resonant layer. If it leaves without any further scattering it has traveled a total distance 
(5t + x + i3t — x = 2(3t inside a resonant layer (see Fig. 3) . It is easy to see that this quantity 
does not change at any scattering is valid for any number of scatterings. Thus, regardless 
of Tres ^ photou speuds a time 2(5thTres/c in the region where it can resonate with electrons. 
This can be compared with scaling oc r for residence time and oc for number of scatter- 
ing for non- resonant scattering. An initial short (5(t)-pulse will, of course, be smeared and 
delayed since depending on scattering history a photon leaves the resonant layer at differ- 
ent spacial coordinates. In the limit of large optical depth spectral smearing is negligible, so 
that all transmitted photons are delayed by j3thTres/c with respect to freely propagating ones, 
(reflected photons are, equivalently, being reflected by a "wall" located at fresi}- ~ Pth))- 

The two most unusual characteristics of resonant cyclotron scattering in a strongly 
inhomogeneous magnetic field are the facts that in the high optical depth limit the transmis- 
sion coefficient is 1 /2 and that the transmitted multiple-scattered photons are preferentially 
up-shifted in energy (for a water-bag distribution they are always up-shifted in energy). Ac- 
cording to (Zhelezniakov 1996), semitransparence of the resonant layer can be explained in 
the following way. In the high optical depth limit, most photons will experience the first 
scatter on electrons moving backwards (along the negative direction x) with (3 ~ —(3th- This 
occurs near the dimensionless coordinate x ~ —(3th- Photons that are re-emitted backwards 
have their frequency increased by 1 + 2(^th- For these photons the resonant layer is located 
closer to the star, shifted by Ax = —2(3th, so that the reflected photon flnds itself outside of 
the resonant layer. Thus, the photons are trapped between two resonant layers. As a result 
their escape probability is 1/2 in each direction (sec Fig. 3). 

To understand the energy up-shift of the transmitted multiple-scattered photons, sup- 
pose that a photon is scattered backward at some position a; on a particle with velocity 
(3i. The energy of the photon increases for /3i < or decreases for /3i > 0. The photon 
can resonate again only with particles whose velocity satisfies (3 > (3i. As a result, in the 
subsequent forward scatter its energy will either decrease by an amount smaller than the 
initial increase (for /3i < 0), or increase by an amount larger than the initial decrease (for 
(3i > 0). 



-14- 



3.4. Validity of approximations 

For non-relativistic plasma, jSth <^ 1, and for low energy photons, e <^ mc^, recoil 
frequency shift is small in each scattering event. Since for resonant scattering number of 
scattering in the limit Tres ^ 1 scales as ~ Tres, the condition for neglect of accumulated 
recoil is ^^T^es <S Pth where e is photon energy. For ~ keV photons and Pth ~ 0.1, this 
limits optical depth to Tres < 10. Since in this case the number of scatterings increases with 
optical depth slower than for non-resonant scattering, this condition on optical depth is less 
stringent. 

For non-relativistic plasma temperatures transition to higher harmonics s are suppressed 
by <^ 1. At high optical depths, tq ^ 1, such transition may still contribute efficiently 
to radiation transfer in magnetic field (similar to the case of energy tranfer in wings of highly 
absorbed lines). The conditions to neglect higher order transition is then tqP^^ <^ 1. 

In case of non-resonant scattering, for high optical depth number density of photons 
increases if compared with the case of no absorption. This leads to increased importance 
of induced processes. For resonant scattering this is relatively small effect since increase 
in photon residence time (and corresponding number density) inside the magnetosphere is 
small, typically a factor 1 -|- (3th ~ 1- 



4. Magnetar spectra due to resonant scattering 

4.1. Magnetospheric plasma in magnetars 

In magnetars, large scale currents flow in the magnetosphere, pushed out of the neutron 
star interior by the slow unwinding of the non-potential magnetic field (Thompson et al. 
2002). The strength of the electrical currents is determined by the speed with which the 
internal magnetic field unwinds. It is expected that the resulting toroidal magnetic field 
Bfj, is limited by the poloidal field Bp, < Bp (for larger fields the configuration is likely 
to be unstable to, e.g. , kink modes). If ~ ^(i>Bp, where < 1 is a pitch angle of 
the magnetic field, the particle density required to carry the current will result in a plasma 
frequency of the order 

ul-^UB--^(t> (40) 
where v < c is a the typical velocity the current carries. The resulting optical depth 

Tres ~ -A0 (41) 
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is of the order of unity, regardless of the mass, frequency and magnetic field (Thompson et al. 
2002). Thus, currents flowing in magnetar magnetospheres may provide large optical depth 
to resonant cyclotron scattering. Electrons then can be heated due to the development of a 
two stream instability operating on short time scales, ~ l/up ^ ^JrJbJBCT- 



4.2. Multiple scattering of radiation in the X-ray band 

Emission from the neutron star's surface may be modified by resonant scattering. As 
radiation from the surface propagates in the magnctosphere it is scattered in a resonant 
layer, some of it is reflected back toward the star and may experience multiple scattering 
between different hemispheres, so that the emergent spectrum is a sum over many reflections. 
At every reflection from a resonant layer radiation is broadened by the random motion 
of resonant electrons, without changing the mean frequency, while the multiply scattered 
transmitted radiation is Compton up-scattered 1 + li^th- If the scattering plasma occupies 
a large fraction of the magnctosphere, the resulting spectrum is a sum of broadened Plank 
spectra (the broader the component, the smaller its relative intensity), Compton up-scattered 
by an optical depth-dependent factor of the order l^lfiih. 

First, we calculate semi-analytically the transmitted fiux for initially Plankian spectra by 
summing over multiple refiections between resonant layers located in different hemispheres. 
In the spirit of our one dimensional approach, we assume that the scattering surfaces are 
plane parallel. If the spectrum of the surface emission is given by ns(a;o) then the observed 
spectrum becomes 

nohs^ J dLLion+{LLi,LLio)ns{coo) + J duJin+{LL! , lji) J n_{LLii,LLio)ns{uJo)duJo -\ (42) 

where the sum is over multiple reflections. Fortunately, the sum quickly converges since 
each consecutive term is at least two times smaller than the previous one (for numerical 
computation we limit the sum to 6 terms, which results in an error of the order of 2% at 
most) . 

Results of the calculations are presented in Fig. 4. First, the spectrum is overall up- 
scattered. The initial spectrum is shifted to higher energies, by as much as 1 -|- 2(3th in the 
limit To » 1 (for smaller optical depths the shift of the emission peak is smaller). As a 
consequence, the real surface temperature is smaller by as much as a factor 1 -|- 2Pth than 
the one inferred from the Compton up-scattered flux. 



Secondly, for optical depth tq ~ 1 the transmitted spectrum is non-Plankian, see Fig. 5. 
Modification of the spectra comes both from increased photon dispersion during each re- 
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flection as well as from summing contributions from several reflections. Note, that the 
non-Plankian form of the transmitted spectrum is mostly prominent at Tq ~ few. In this 
case the reflected flux becomes comparable to the initial flux. On the other hand, in the 
limit of large optical depth tq ^ 1, each reflection does not lead to an increase in the flux 
dispersion, so that the transmitted flux remains mostly Plankian, but with the peak shifted 
to larger energies by 1 + 2(3th- 

The asymptotic spectrum at frequencies much higher than the initial frequency can be 
estimated if we assume that each reflection acts as a Gaussian filter with a typical dispersion 
c ~ A/i^o < ^0 and that the reflection probability is < p < 1/2. Neglecting the small 
frequency shift during the final transmission, the resulting flux F is a sum over all reflections 

F oc V ^ . (43) 

For n ^ 1 we can substitute the sum by an integration and evaluate the integral by the 
steepest decent method at n ~ uj / {^/2AIl{TJp)a) , 



1 



F oc (44) 

Qualitatively, it takes on average n ~ uj/{f3th^Q) scatterings to reach a frequency u ^ ujq, 
but the intensity of radiation is decreased by p""*"^. Thus, scattering by weakly relativistic 
warm plasma with (3th < 1 produces a non-thermal spectrum with an exponential cut-off at 
uj ~ kBT/{h(5th)- 



4.3. Spectral Fitting 

Motivated by the fact that resonant Compton scattering qualitatively reproduces mag- 
netar spectra, we performed a preliminary analysis to verify that it can also do so quantita- 
tively. We decided to flt our model to a 28.8 ks long Chandra X-ray observatory observation 
of the anomalous X-ray pulsar (AXP) IE 1048.1-5937. The data were taken with the ACTS 
detectors (S3 chip) in CC mode. In the literature IE 1048.1-5937's spectrum has been charac- 
terized by a two-component model consisting of a blackbody of temperature ksT ~ 0.6 keV 
and a power-law tail with spectral index F ~ 3. We chose this source because its spectrum 
is intermittent between the anomalous X-ray pulsars (AXPs) and the soft-gamma repeaters 
(SGRs). An analysis of all magnetar spectra will be done in the future. 

Given that our model is semi- analytical we created a table of intensities for different 
values of Pth, Tres ^ud ksT. We developed software to convert this grid into a format 
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that could be read by the X-ray spectral fitting program XSPEC^. Using XSPEC the model 
was folded through the instrumental response and multiplied by a photoelectric absorption 
model. Only counts between 0.5-1.67 and 1.91-6.7 keV were included in the fits, leaving 281 
channels for spectral fitting. The 1.67-1.91 keV region of the spectrum was ignored because 
of an instrumental feature in the spectrum. For more details on the analysis of this Chandra 
observation of IE 1048.1-5937 see (Gavriil et al. 2005). In our fits jSth and Tres were held 
fixed at predetermined values but we allowed the column density Nh, the temperature ksT 
and the normalization to vary freely. We repeated the procedure for different values of (5th 
and Tres until we minimized a statistic. Even with our course grid in (5 and r^es we were 
able to successfully reproduce the observed spectrum. The fit had a reduced x^^f = 1.05 for 
278 degrees of freedom (dof). We compared our model to a blackbody plus power law model 
and a double blackbody model. A summary of our spectral fits is given in Table 1. 

For comparison we plot the best-fit blackbody plus power law model and the best-fit 
resonant Compton scattering model for this observation in Figure A. It is interesting to note 
that our model fairs just as well as the "canonical" magnctar spectral model of a blackbody 
plus power-law model (see Table 1). Furthermore, it is interesting to note that all three 
models listed in Table 1 have the same number of parameters and fit the data well, however, 
the resonant Compton scattering model has an advantage: it is a single (rather than a 
double) component model which has a clear physical interpretation within the framework of 
the magnetar model. 



4.4. Seculeir spectral vairiations in magnetars 

Most magnetar candidates (AXPs and SGRs) show evidence of secular spectral variabil- 
ity {e.g. Marsden & White 2001; Woods & Thompson 2004). In the magnetar model this 
is readily explained by a variable current fiowing in the magnetosphere producing variable 
optical depth to resonant scattering (Thompson et al. 2002). ^ Spectral variations in the 
persistent emission of magnetar candidates have been prominent during active bursting pe- 
riods. For example, significant spectral evolution was observed during an outburst of AXP 
IE 2259.1-1-586 involving ~80 bursts within a 3 hr observations (Kaspi et al. 2003). More 



^http: / /xspec. gsfc.nasa.gov 

^Isolated neutron star RX J0720. 4-3125 also shows spectral variations of its black-body component on 
time scale of years Paerels et al. (2001); de Vries et al. (2004). Resonant scattering may offer a possible 
explanation, but an upper limit on the magnetic field < 3 x 10^^ G Kaplan et al. (2002) seem to exclude 
the magnetar hypothesis for it. 
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recently a similar, yet much more subtle change was reported by Rea et al. (2005) in AXP 
IRXS 1708-4009. They find that the source's spectrum has been softening, specifically they 
find that the photon index increased from V ~ 2.4 to ~2.8. In both cases the correlation 
between the luminosity and the hardness of the non-thermal component of the source's spec- 
trum is consistent with the large-scale current bearing magnetosphere model of Thompson 
et al. (2002). 

5. Other applications 

5.1. Smearing of spectral lines 

Atmospheres of neutrons stars were expected to produce spectra abundant with spec- 
tral features [e.g. Zavlin & Pavlov 2002). Though spectral features were observed in some 
pulsars, their absence in some prominent sources hke RX J1856. 5-3754 (Drake et al. 2002), 
poses a challenge to atmospheric models, (see, e.g. , van Kerkwijk 2003, for review). Reso- 
nant scattering may provide a possible explanation. Let us calculate the effects of cyclotron 
scattering in a magnetosphere filled with plasma on a narrow spectral feature superimposed 
on a thermal spectrum. The initial spectrum is taken to be a thermal plus a narrow emission 
line located at the energy eum = ^ksT with a width Aeune = O-lksT and total intensity 
1/100 of the blackbody (Fig. 6). The emission passes through the magnetosphere with optical 
depth tq — 1. After multiple reflections the transmitted intensity is virtually indistinguish- 
able from the one obtained from just a thermal spectrum without an emission line. Thus, 
scattering in the magnetosphere may provide an explanation for the lack of spectral features 
in neutron star. 



5.2. Optical emission 

The character of resonant scattering changes for UV and longer wavelengths, when the 
wave polarization at resonance becomes circular Eq. (12). Consider, for simplicity, the scat- 
tering of unpolarized light in a monopole magnetic field in a one-dimensional approximation 
in the case when normal modes are circularly polarized. For unpolarized light propagating 
through an electron-ion plasma, half of the intensity (with circular polarization such that 
electric field of the waves rotates in the same direction as the positive charges) is transfered 
without scattering, while the other half (with polarization such that the electric field rotates 
with the electrons) is scattered with the transfer functions calculated in Section 3. In the 
limit Tres ^ 1 half of the scattered radiation is re-emitted forward, half backward, so that 
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the transmission coefficient is 3/4 instead of 1/2 for the case of hnearly polarized waves. 

In addition, two possibihties exist for multiple scattering of a reflected photon from 
the opposite hemisphere. If a photon propagates mostly outside the polarization limiting 
radius, its polarization does not change along the ray, while for propagation inside the 
polarization limiting radius its polarization evolves adiabatically. These two possibilities 
result in qualitatively different scattering regimes: if a photon keeps its polarization, then 
in the simple one-dimensional picture in the opposite hemisphere the electric field of the 
wave rotates in the opposite direction to that of the local electron so that photon will not 
experience any more scatterings. On the other hand, if polarization evolves adiabatically, 
by the time a photon would reach a resonance in the opposite hemisphere it will become 
virtually unpolarized and thus can be scattered. 

Location of polarization limiting radius rip can be estimated from the condition Anrip ~ 
A where A is wavelength and An is the difference between refractive indexes of the two normal 
modes. In the relevant regime An = ^{B / BqedY (Heyl & Hernquist 1997), which gives 



where a fin is fine structure constant. This is typically larger than resonance radius, so that 
polarization varies adiabatically. 



The importance of cyclotron scattering in the magnetosphere on the soft X-ray proper- 
ties of young pulsars has been first pointed out by Wang et al. (1998), but see also Ruderman 
(2003). It was suggested that abundant pair creation may occur on closed field lines which 
would produce a large optical depth to cyclotron scattering. Unfortunately, the model of 
the cyclotron blanket proposed in the above references is incorrect since it was assumed 
that in the high optical depth regime most of the emission is refiected back to the star. In 
fact, resonant up-scattering leads to an increase of the radiation temperature by a factor 
~ 2Pth ~ 10 — 20% at most, contrary to what was proposed by Ruderman (2003). 

On the other hand, estimates of the effective emitting area based on the observed 
luminosity and the effective temperature do change due to resonant scattering. Radii 
of neutron stars are conventionally found from the luminosity and observed temperature, 
Rbb ~ {Lx / ^ncTsBT^y^'^ , where Lx is the observed X-ray luminosity, (Jsb is Stefan-Boltzmann 
constant and T is the temperature inferred from the spectral fit. This is called blackbody 
radius. Resonant scatering changes estimates of both temperature and luminosity and thus 




(45) 



5.3. Young pulsars 
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changes estimate of neutron star radius. If effects of resonant Comptonization are important, 
then the surface temperature is lower than the observed one: Tg ~ Tx/(1 + 2aPth), (where 
< a < 1 is a coefficient dependent on the resonant optical depth). On the other hand, 
the observed luminosity has contributions both from the thermal surface emission and from 
the energy of the scattering electrons. In case the number of photons is conserved the total 
luminosity is Lx ~ (1 + 2a(3th)Lx,s- Thus, the real radius is larger than blackbody radius: 

Rns ~ iLx,s/^7rasBT^y^'' = Rbbil + ^a^thf'' (46) 

This may have important implications for attempts to determine equation of state of neutron 
stars. 



5.4. Accreting neutron stars 

For sources accreting through the magnetosphere at a residual rate M, the optical depth 
at energy e is 

neM I M \ .y^ / e \ i7/6 



^/rresGMrUpBres \ lO-i^Mo/yr 



with Tres givcu by Eq. (1). This accretion rate is larger than the Bondi rate for a neutron star 
moving at ~ 100 km/sec through the ISM, Mb ~ lO^^''M0/yr, so that isolated neutron stars 
which are not high energy pulsars or magnetars cannot support dense, hot magnetospheres 
(cold plasma may be supported by radiation pressure Dermer & Sturner (1991)). On the 
other hand, the accretion rate (Eq. 47) is considerably smaller than is typically reached in X- 
ray binaries. It is expected that most of the accreting material is either channeled along the 
polar magnetic field lines or is flung out by the propeller effect, but even if a small amount 
of accreting material gets into the bulk of the neutron star's magnetosphere {e.g. due to 
development of resistive instabilities in the boundary layer) the resulting resonant cyclotron 
optical depth can be considerable. 

In addition, even in the accretion column, the effects of magnetic field inhomogeneities 
are non-negligible since the size of the accretion column H, which is typically tens of meters, 
may become comparable to the width of the resonance layer, which for temperatures of ~ 5 
keV {I3th = 0.01) and dipole magnetic field Lb ~ f^Ns/^ is Sr = PthLB ~ 30 meters. Thus, 
effects of magnetic field inhomogeneities should introduce corrections of the order of unity 
to resonant radiative transfer in the accretion column. 



- 21 - 



6. Discussion 

An astrophysically small amount of matter, heated by some external source of energy, 
may considerably modify the emergent spectra through resonant cyclotron scattering in the 
magnetospheres of neutron stars and thus impede the effort to understand the structure 
and composition of neutron stars and their equation of state. The way emission is modified 
is qualitatively different both from the case of non-resonant Thomson scattering as well as 
resonant scattering in homogeneous magnetic field . 

During resonant Comptonization the number of photons is conserved, while their energy 
increases by a factor of ~ 1 + 2Pth- This photon energy increase comes at the expense of the 
electron energy. Since the minimal charge densities required to produce an optical depth of 
the order of unity are small (Eq. 9), leptons are efficiently cooled by Comptonization, so that 
to be in steady state they should be heated by some energy source. The power required is 
~ PthLx where Lx is the observed X-ray luminosity. Depending on the nature of an object 
heating can be achieved by different sources. In young pulsars this may be achieved by 
pair production by high energy 7-ray emission (Wang et al. 1998). In magnetars, charges 
are pulled from the surface by inductive electric fields generated due to slow crustal motion 
(Thompson et al. 2002). In X-ray binaries residual accretion through the magnetosphere 
will heat the plasma due to the release of gravitational energy. 

Resonant Comptonization by mildly relativistic plasma generically results in soft non- 
thermal spectra, similar to what is observed in anomalous X-ray pulsars at intermediate 
energies, < 10 keV. Non-thermal effects are most prominent for mildly relativistic plasma 
Pth < 1 at intermediate optical depths, r^es ~ 1- On the other hand, at large optical depths, 
Tres ^ 1 transmitted spectrum is nearly thermal, but shifted to higher frequencies. 

An important advantage of the model is that in our case photons are just moved to 
higher energies, so that one cannot simply extrapolate the data from the X-rays down to 
longer wavelengths. In the conventional blackbody plus power law fitting procedure one also 
needs to assume a low energy cut-off to the power law, (c/. Halpern & Gotthelf 2005). Such 
a cutoff is not required in our model. 

As a principal application of the model, we have fitted the non-thermal spectrum of 
AXP IE 1048.1-5937 and found a quantitatively good agreement. Besides the production of 
non-thermal spectra in magnetars, resonant scattering may change estimates of the surface 
temperature while conserving the overall Plankian form of the spectrum (in case of high 
optical depths), and may explain the lack of surface spectral lines. Variations of magneto- 
spheric plasma parameters may induce spectral changes on time scales much shorter than 
the evolutionary time scale of neutron stars. 
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In this paper we addressed the principal issues related to resonant cyclotron transfer, 
employing a simple one- dimensional model of resonant radiation transfer in the limit of non- 
relativistic temperatures of scattering particles with a water-bag distribution. More involved 
three-dimensional Monte- Carlo simulations are in progress. They will extend the current ap- 
proach to polarized surface emission. The magnctospheric plasma may be relativistically hot, 
(3th ~ 1, in which case the diffusion of photons both in frequency and in real space is nonlo- 
cal and, in addition, scattering on high harmonics becomes important. Also, possible bulk 
motion of the scattering plasma may lead to the additional effect of bulk Comptonization. 
Resonant scattering by relativistic electrons has been invoked by Thompson & Beloborodov 
(2004) to explain the high energy emission from magnetars. Finally, the fact that surface 
radiation is not blackbody {e.g. Zavhn & Pavlov 2002) should also be incorporated. We 
plan to address some of these issues in a subsequent publication. 

We would like to thank Victoria Kaspi and Christopher Thompson for numerous dis- 
cussion and insightful comments. We arc grateful to Peter M. Woods for his help reducing 
the Chandra observation of IE 1048.1-5937. We also would like to thank W. Ho, K. Mori, 
M. Ruderman, M. van Kerkwijk for comments and discussions. 
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A. One-dimensional Monte-Carlo simulation 

Analytical results in Section 3 have been obtained in the limit /3 — > 0, neglecting the 
kinematic enhancement of the scattering rate by a factor 1 — P due to the motion of a 
particle. To test the validity of this approximation we perform one-dimensional Monte Carlo 
simulations of the transfer of resonant cyclotron radiation. In the limit /3 <S 1 we can expand 
frequencies about the initial frequency {cu — ujq)/ujq <S 1 and assume ujb — ojq{1 — x) where 
X is a dimensionless coordinate (measured in terms of L^). Under these circumstances the 
resonant velocity is /? = ±(a; + uo — loq) and at each scattering the new frequency becomes 
Co = 2ujq{1 — x) — uj. At each step the probability of scattering backwards is multiplied by a 
kinematic factor 1 — j3. Results for the transmission of a spectral line arc present in Fig. 8. 
At intermediate optical depths back scattered photons have on average higher frequencies. 
This is due to the fact that photons are more likely to encounter backward propagating 
electrons, whose cross-section is enhanced by 1 -\- (3. On the other hand, at high optical 
depth deviations from analytical results are small. In this limit photons experience multiple 
scatterings inside the layer so that the relative velocity of the electrons between each scatter 
is small and the kinematic enhancement is not important. 



This preprint was prepared with the A AS IM^jX macros v5.2. 
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Table 1. Phase-averaged spectral fits of IE 1048.1-5937 for various models. 



Parameter Value* 



Resonant Compton Scattering Model 



Nh (10^2 cm-2) 




0.76 ±0.16 


ksT (keV) 




0.39 ±0.05 


P 




0.3 (fixed) 


'^res 




3 (fixed) 


Flux*^ (10-^2 gj-gs cm-2 s-i) 




7.61 


Unabs Flux^ (lO"^^ ergs cm"^ s" 




8.64 


xVdof 




293/278 


Blackbody + Power Law 


Nh (1022 cm-2) 




1.18 ±0.04 


keT (keV) 




0.585 ±0.012 


r 




3.08 ±0.10 


FIUX*^ (10~^2 g2;gg cm~2 s~^) 




7.91 


Unabs Flux*^ (10~^2 ^^gs cm~2 s" 




9.56 


xVdof 




274/276 


Blackbody + Blackbody 


Nh (1022 cm-2) 




0.76 ±0.03 


ksTi (keV) 




0.458 ± 0.030 


keT^ (keV) 




0.843 ± 0.058 


Flux''(10~^2 ergs cm-2 g-i^ 




7.61 


Unabs Flux^ (10"^2 ergs cm-2 g- 




8.62 


xVdof 




284/276 
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^The quoted uncertainties reflect the la error for a reduced of unity. 
''Flux in the 2-10 keV band. 
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-3Po -Po Po ^ 

Fig. 1. — Region of resonant interaction in {x, (3) and (,^, rj) coordinates (77 = f3—x, = x+f3) 
for a water-bag distribution. Monochromatic light with frequency ujq propagates through 
the resonance layer centered at usix = 0) = ujq and is subsequently scattered by particles 
occupying the lightly shaded area \P\ < (3th- Resonance transfer occurs in the darkly shaded 
triangle delimited by the lines 77 = 0, ^ = —'^Pth, V — "^Pth — 
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Fig. 2. — Transmitted (solid lines) and reflected (dashed lines) fluxes for a water-bag dis- 
tribution of particles with f3tfi = 0.1 and different optical depths Tq. For small Tq < 1 the 
reflected photons have a homogeneous distribution since the intensity of the initial beam re- 
mains approximately constant. The reflected intensity is oc tq. The transmitted flux shows 
an initial 5- function and a double scattering tail oc r^, the different panels show that as the 
optical depth increases the height of the initial (5-function peak decreases, while the intensity 
of reflected beam increases. In the limit r^es — oo both the reflected and the transmitted 
intensities are equal (one half of the initial intensity) . For larger optical depth the dispersion 
of both the reflected and transmitted intensities decrease . In this limit the transmitted 
photons are Doppler up-shifted on average by Acu/lu = H- 2Pth, while the reflected photons 
have on average a; ~ a;o. 
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Fig. 3. — Illustration why photon always spends dimensionless time 2(3th in the resonance 
region {2l3thLB in regular units, where is the magnetic field inhomogenuity scale) and why 
resonant layer is semi-opaque in the limit of high optical depth. Initially, resonant surface 
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Fig. 4. — Modification of an initially Plankian spectrum (dashed line) by multiple cyclotron 
scattering; Pth — 0.3 for different values of Tres- 
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Fig. 5. — A fit to the transmitted spectrum (solid line) for optical depth tq = 3 and f3th = 0.3. 
Dashed line is the initial Plank spectrum with Tq = 1, Dotted lines are a Plank fit to the 
transmitted spectrum at the peak of the emission, T — 1.25 and a residue. The residue is 
strongly non-thermal. 
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Fig. 6. — Effects of resonant scattering on a line in a spectrum. Initial spectrum (dashed line) 
is Plankian plus an emission line located at hwune = 5k bT, with a width Auum = 0.1 /c^r 
and total intensity 1/100 of the black body. The transmitted intensity (solid line) is barely 
above the transmitted intensity for the no line case, Fig. 4, (dotted line). The difference 
between the two (long dashed line) has a complicated profile, but it's intensity is only ~ 1%. 
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Fig. 7. — X-ray Spectrum of the Anomalous X-ray Pulsar IE 1048.1-5937 as observed by the 
Chandra observatory. Left: The data fit with our resonant Compton scattering model. See 
Table 1 for best-fit parameters. Right: Same spectra but fit to a blackbody plus power law. 
see Table 1 for the best-fit parameters. The dashed and dotted hnes represent the individual 
model components. See text for details on the analysis 
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Fig. 8. — Same as Fig. 2 reproduced using one-dimensional Monte-Carlo simulations with 
10^ photons and taking the 1 — /9 enhancement of the cross-section into account, Pth — 0.1. 
At intermediate optical depths the reflected flux is skewed towards higher frequencies due 
to the kinematic enhancement of the cross-section. These simulations are done using a 
one-dimension approximation for constant Lb- 



